平 谱 射电 类 星体 在 着 康 普 顿 机 制 中 软 光子 来 源 的 研究 * 
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摘要 : 为 了 验证 逆 康 普 顿 机 制 中 的 主导 过 程 ， 以 及 主导 过 程 中 软 光 子 的 起 源 。 我 们 通 
过 文献 收集 了 75 个 平 谱 射电 类 星体 多 波段 的 ( 准 ) 同时 性 数据 ， 并 利用 对 数 抛物 线 函 数 拟 合 
了 频谱 能 量 分 布 , 分 析 同 步 峰 流量 与 逆 康 普 顿 恬 之 比 和 多 普 勒 因子 的 关系 , 运用 自 举 法 计算 
了 样本 的 种 子 因子 分 布 用 以 调查 软 光 子 的 起 源 。 分 析 结 果 如 下 : (1) 、 在 轻 子 模型 框架 下 ， 
整个 样本 与 子 样本 log(UBp(vE,)Ec/Uext(vVE,)syn) 与 og6 的 显著 相关 性 表明 平 谱 射 电 类 星体 的 
高 能 辐射 优先 由 外 康 普 顿 过 程 主导 。 (2) 、 对 于 相对 低 同 步 峰 峰 频 的 种 群 ， 观 测 的 种 子 因子 
分 布 与 理论 上 和 尘埃 环 的 种 子 因子 范围 相交 表明 其 逆 康 普 顿 过 程 中 的 软 光子 很 可 能 由 尘埃 环 
主导 , 而 对 于 相对 高 同步 峰 峰 频 的 种 群 ， 观测 的 种 子 因 子 超过 了 宽 线 区 或 尘埃 环 种 子 因子 的 
范围 ， 因 此 我 们 还 无 法 确定 其 软 光子 的 来 源 。 
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0 引 言 


WEAR PK (Blazar) 是 观测 者 视线 沿 相 对 论 性 喷 流 轴 向 方向 观测 到 的 射电 噪 活动 星系 核 品 
(Radio-Loud Active Galactic Nuclei) 的 一 个 子 类 ， 包括 可 以 观测 到 强 发 射线 的 平 谱 射 电 类 
星体 (Flat-Spectrum Radio Quasar, FSRQ ) 和 观测 到 弱 或 者 无 发 射线 的 蝎 虎 天 体 (BL Lacerate, 
BL lac) 两 个 子 类 口 。 观 测 上 ， 耀 变 体 的 宽 波 段 辐射 是 从 射电 到 伽 马 射线 波段 的 光谱 能 量 分 
fü (Spectral Energy Distribution, SED) ， 在 logv -logz 玉 图 上 呈现 出 典型 的 双 峰 结构 。 前 
人 对 这 个 双 峰 结构 进行 了 大 量 的 研究 ， 在 目前 应 用 最 为 广泛 的 轻 子 模型 框架 下 ， 低 能 峰 被 
认为 是 相对 论 电 子 的 同步 辐射 (Synchrotron Radiation, syn) ， 高 能 峰 被 认为 是 相对 论 电 子 
与 软 光 子 作用 而 产生 的 逆 康 普 顿 散射 名 (Inverse Compton cattering, IC) 。 如 果 软 光子 来 自 
同步 辐射 ， 称 为 同步 辐射 自 康 普 顿 模型 中 (Synchrotron Self-Compton, SSC) ; 如 果 软 光子 
来 自 喷 流 外 部 ， 则 称 为 外 部 康 普 顿 模型 (External Compton, ЕС) 。 在 外 康 普 顿 模型 的 情况 
To ОКР n] ВЕЖУ FIKIRA (Accretion Disc, AD), 3: ££ x^ (Broad-Line Region, ВІВ) 
BYP IA" (Dust Torus, DT) 。 到 目前 为 止 ， 逆 康 普 顿 散 射 的 软 光 子 来 源 仍 然 无 法 完全 确 
ce Foe, 


È 


“基金 项 目 : 国家 自然 科学 天 文联 合 基金 〈U1431123，10973034) 资助 . 
WS AHA: 2021-01-18; ”修订 日 期 : 
作者 简介 : 刘 小 脑 ， 


通讯 简介 : 谢 照 华 ， 


硕士 . 研究 方向 : 天 体 物 理 . Email: boyliuxiaopeng@ foxmail.com 


dB 8 


， 教 授 . 研究 方向 : 理论 物理 和 天 体 物 理 . Email: zhxie007@126.com 


Chinay n "M" НА | 
Chinax ive КЕЛНУ 


一 些 研 究 者 认为 耀 变 体 的 短 时 标 光 变现 象 意味 着 其 辐射 产生 位 置 在 BLR 以 内 "在 
标准 的 单 区 轻 子 模型 框架 下 , 无 法 解释 FSRQ ЈЕ ЗС 279! 3C 454.3"? DL 4C 21.35 5 
SUA ЕВЕ] Cc 100GeVO 现象 。 利 用 其 长 基线 干涉 技术 5 CVery-Long Baseline 
Interferometry, VLBI) 对 射电 辐射 进行 观测 , 研究 者 们 发 现在 射电 波段 有 一 些 明亮 的 区 域 ， 
这 些 区 域 被 称 为 射电 核 ， 其 位 置 远离 中 心 黑洞 ， 其 至 在 DT 之 外 中 。 文 [16]Tan 收集 了 60 
个 有 多 波段 准 同时 性 数据 的 FSRQs, 在 单 区 模型 框架 下 , IC 机 制 中 同时 考虑 了 BLR 和 DT 
贡献 的 种 子 光子 ， 他 们 发 现 FSRQs 的 高 能 辐射 来 源 于 BLR 之 外 DT 之 内 。 文 [17]Harvey 
收集 了 62 个 有 多 波段 准 同时 性 数据 的 FSRQs， 在 轻 子 模型 的 框架 下 ， 引 入 一 个 参量 
于 观测 到 的 多 波段 SEDs 去 计算 双 峰 结构 的 峰 频 频率 和 通 量 ， 来 确定 blazar 中 IC 机 制 下 
EC 过 程 中 主要 的 软 光 子 场 。 他 们 把 这 个 参量 称 为 种 子 因 子 〈Seed Factor, SF) 。 结 合 运用 
自 举 法 (Bootstrapping) 发 现 观测 的 种 子 因子 与 期 望 上 来 自 DT 的 种 子 因子 在 统计 上 显著 一 
致 ， 与 期 望 上 的 BLR 不 一 致 。 他 们 研究 结果 表明 IC 机 制 中 的 软 光 子 场 主要 由 DT 贡献 。 

综 上 所 述 ， 目 前 仍然 无 法 确定 IC 辐射 的 软 光 子 起 源 ， 也 无 法 确定 辐射 产生 的 位 置 。 在 
本 文中 ， 我 们 使 用 包含 同时 性 的 大 样本 从 统计 视角 分 析 IC 机 制 中 软 光 子 的 起 源 。 


== 
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1 样 本 


本 文 从 Feimi 4-LAC 活动 星系 核 目 录 收 集 了 69 个 有 多 波段 准 同时 性 数据 的 平 谱 射电 
类 星体 中 ， 另 外 在 文献 [19] 中 收集 了 6 个 有 多 波段 准 同 时 性 数据 的 平 谱 射 电 类 星体 样本 ， 
E75 个 平 谱 射 电 类 星体 。 我 们 收集 的 样本 都 满足 下 列 条 件 : 有 确定 的 红 移 ; 不 考虑 射电 波 
段 的 数据 ， 履 盖 紫 外 到 伽 马 射 线 波 段 ， 各 个 波段 相交 的 观测 时 间 小 于 2 个 月 ， 可 以 构建 双 
峰 结 构 的 SEDs;， 有 可 以 查 到 的 体 洛 伦 效 因子 名 或 多 普 勒 因子 59。 在 本 文中 ， 我 们 假设 体 
洛 伦 效 因子 与 多 普 勤 因子 近似 相等 e， 使 用 的 磁场 СВ) 参数 由 文 [16] 中 的 模型 拟 合 得 到 。 
© 我 们 的 样本 有 一 个 相对 高 峰值 频率 的 源 (PKS 21414175, Товугуп = 15.08 Hz) . 

分 别 使 用 ASDC SED Builder Ж Н Кети ВЕЖ SEDL”, Python 拟 合 提取 自 文 [19] 
A RA BU SEDs. ФИТАН RERI Ns HER AS SED 的 低能 峰 和 高 能 峰 ， 二 次 函数 
RADA”, 


log(vFv) = c(logv)? + b(logv) + а (1) 
其 中 ，a、b 和 c 为 常数 。 我 们 通过 拟 合 观测 数据 得 到 的 峰值 频率 和 整个 峰 的 积分 光度 都 是 
在 观测 者 坐标 系 下 的 , 为 了 去 除 红 移 的 影响 , 我 们 需要 将 这 两 个 量 转 到 AGN 静止 参考 系 ， 
转换 公式 如 下 : 


Vsynic = (1+ Хуга с (2) 
(VLy) sync = AnD? (УА,) увис (3) 
X OD 中 , VSic 是 观测 坐标 系 下 的 同步 峰 或 逆 康 普 顿 峰 峰 值 峰 频 ，z 为 红 移 。 式 (3) 中 ， 


? https://tools.ssdc.asi.it/index.jsp 


Di AGEIA, (VE, )synrc EMMA A hs ARP ES [8] 25 ЈА ВЈ EHE DUAE DP ВА E о 

表 1 中 展示 了 我 们 收集 的 75 个 样本 , 第 (1) 栏 , 源 的 名 称 ; 第 (2) 栏 , 红 移 ; 第 (3) 
栏 ，( 准 ) 同时 性 时 间 〔( 文 [19] 中 的 数据 未 注 明 时 间 〉; 第 CAO 栏 ， 取 对 数 后 的 同步 峰 峰 
值 频率 (以 Hz 为 单位 ) ; 第 (5) 栏 ， 取 对 数 后 的 同步 峰 流量 (以 erg cm? s 为 单位 ); 
第 (6) 栏 ， 取 对 数 后 的 同步 峰 光 度 〔 以 erg s 一 为 单位 ) 。 第 〈7) 栏 ， 取 对 数 后 的 IC WE 
峰值 频率 〈 以 Hz 为 单位 ) ; 第 (8) 栏 ， 取 对 数 后 的 逆 康 普 顿 峰 流 量 〈 以 erg cm"? s 1 为 
单位 ) ; 第 (9) 栏 ， 取 对 数 的 逆 康 普 顿 峰 光 度 〈 以 erg s 一 为 单位 ) ; 第 (10) №, 磁场 
能 量 密度 (以 erg 5 为 单位 ) ; 第 (11) 栏 , 体 洛 伦 兹 因子 (多 普 勒 因子 以 上 角 标 a 注 明 ); 
第 (12) 栏 ， 观 测 的 种 子 因 子 。 


Ц 


#1 样本 


Table.1 The sample 


Source name Z Time logvsyn log(vF,)syn log(vLy)syn logvic  log(vE)ic  log(vLy)ic Ug r SF 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (1) (12) 
PKS 0537-286 3.104 2010/3/7-12 12.10 -11471 47.453 20.472 -9.541 49.383 0.05218 12.73 5.599 
АС +04.42 0.966 2010-6-17 12.41 -10.738 46.943 21.749 -9.860 47.821 0.05361 12.11 4.020 
TXS 0653-033 0.634 2010-3-31 2.58 -10.991 47.350 22.854 -10.201 47.009 0.07689 22.49 3.119 
PKS 0605-085 0.872 2010-6-7 12.62 -10.914 46.635 22.989 -10.383 47.167 0.10023 14.94 2.904 
3C 279 0.536 2010-1-15 2.62 -9.764 48.565 22.141 -9.173 47.861 0.01743 12.82 3.694 
PKS 2325+093 1.841 2010/6/18-29 12.64 -11.133 47.219 21.546 -9.780 48.572 0.17865 26.84 4710 
PKS 0215+015 1.715 2010/1/29-2010/2/2 2.67 -10.419 47.878 21.895 -9.947 48.350 0.13171 25.67 3.908 
PKS B1921-293 0.353 2010-9-30 12.73 -9.742 50.088 21.151 -9.764 46.840 0.03922 11.33 4.467 
S5 1357+76 1.585 2010/5/24-2010/5/25 2.73 -11.318 48.842 22.267 -10.556 47.635 0.08655 14.34 3.763 
PKS 1244-255 0.635 2010-6-25 12.74 -10.472 46.761 22.164 -9.462 47.771 0.03841 10.7 4.007 
PKS 1730-13 0.902 2010-3-14 2.75 -10.218 47.389 22.946 -9.190 48.418 0.02501 13.78 3.261 
3C 446 1.404 2010/5/22-27 12.75 -10.383 47.678 21.818 -10.028 48.033 0.04872 12.00 4.048 
3C 454.3 0.859 2009/12/4-6 2.76 -9.234 48.299 22.029 -8.210 49.324 0.19190 16.95 4.173 
3C 273 0.158 2010/6/7-23 12.76 -9.478 46.334 20.247 -8.575 47.239 0.04930 20.87 5.947 
53 0827424 0.942 2010/4/18-19 2.77 -10.595 47.512 20.788 -11.064 46.569 0.09886 13.71 4.668 
NRAO 512 1.66 2010-8-7 12.78 -10.879 47.383 22.800 -9.947 48.315 0.07356 14.87 3.370 
B2 0743425 2.979 2010-10-15 2.78 -11.595 47.451 21.697 -10.192 48.668 0.04276 13.48 4.722 
TXS 0800+618 3.033 2010/4/3-4 12.79 -11.479 47.400 21.626 -10.339 48.539 0.48790 15.11 4.691 
PKS 0454-234 1.003 2010/2/24-25 2.80 -10.417 47.304 22.122 -10.093 47.628 0.03834 8.66 3.764 
4C 38.41 1.814 2010-3-7 12.82 -10.377 49.673 22.082 -9.194 49.142 0.06195 15.33 4.302 
4C+47.44 0:735 2010/7/30-2010/8/1 2.85 -10.764 47.411 21.808 -10.093 47.274 0.03971 9.64 4.280 
ОР 313 0.997 2009/12/12-21 12.86 -10.530 47.163 22.323 -9.900 47.795 0.19773 31.22 3.795 
ОК 290 0.707 2010/5/7-2010/6/15 2.88 -10.721 47.734 21.545 -10.228 47.098 0.12912 14.44 4.497 
В2 1324+22 1.4 2010-3-30 12.88 -10.959 47.100 21.911 -10.080 47.978 0.06439 12.06 4.346 
PKS 2209+236 1.125 2009/4/15-21 2.92 -10.886 46.937 23.142 -10.488 47.335 0.01172 9.30 2.817 
0Ј535 1.418 2010/3/28-2010/4/8 12.93 -11.219 46.853 21.506 -10.357 47.715 0.05034 12.16 4.793 
PKS 1335-127 0.539 2010/1/18-26 2.94 -10.237 47.553 21.367 -10.035 47.004 0.12695 22.62 4.567 


TXS 2331+073 

OI 280 

B3 0908+41 

PKS 1144-379 

54 0917+44 

OG 505 

TXS 0404+075 
PKS 1954-388 

PKS 1546+027 
BZQ J1514+4450 
B2 2319+31 

PMN J2135-5006 
PKS B1434+235 
PKS 2345-16 

PKS 0436-129 

B3 0650+453 

PKS 2255-282 

PKS 0541-834 

B2 2155431 

B2 1520+31 

4C +28.07 

PKS 1406-076 

PKS 2155-83 

PKS 1502+106 

B2 1846+32A 

B3 1708+433 

PMN 2345-1555 
PMN J1959-4246 
PKS 0420-01 

PKS 1203-26 

PKS 0906+01 

PKS 1510-089 
MG2 J110606+2812 
4С+29.45 

MG2 Ј071354+1934 
РММ Ј2145-3357 
PKS 0601-70 

ОС 457 

PKS B1035-281 
PKS 0347-211 
CRATES J1208+5441 
PMN J1347-3750 
MG2 J153938+2744 
CLASS J1333+5057 


0.401 
0.889 
2.563 
1.048 
2.19 
1.254 
1.133 
0.63 
0.414 
0.57 
1.489 
2.181 
1.548 
0.576 
1.285 
0.928 
0.926 
0.774 
1.488 
1.487 
1.213 
1.494 
1.865 
1.839 
0.798 


1.027 
0.621 
2.178 
0.916 
0.789 
1.024 
0.36 
0.843 
0.724 
0.54 
1.361 
2.409 
0.859 
1.066 
2.944 
1.344 
1.3 
2.191 
1.362 


2009/12/20-21 
210/4/1-12 
2010-2-25 
2010-6-24 
2009/10/28-30 
2010-4-25 
2010/2/14-15 
2010/4/9-14 
2010/2/13-2010/2/20 
2010-4-6 
2009-5-20 
2010/4/22-/2010/5/5 
2010-6-14 
2009/12/4-5 
2010-7-1 
2010/3/22-24 
2010/5/20-26 
2009/11/8-11 
2009/7/8-12 
2010-2-5 
2010-5-23 
2010/7/5-17 
2010-7-29 
2010/10/6-19 
2009/12/1-2 
2010/4/1-30 
2010/4/5-14 
2009-8-27 
2010/11/1-2010/12/31 
2010/4/18-2010/5/5 
2010/2/07-2010/3/07 
2010/5/28-2010/6/11 
2009/12/15-17 
2009/9/22-24 
2008/12/1-30 


2010/1/22-23 

2010/5/16-21 
2010-9-20 
2010-3-17 

2010/3/12-17 


2.94 
12.95 
2.97 
12.97 
2.97 
12.98 
2.99 
13.02 
3.04 
13.05 
3.05 
13.05 
3.06 
13.07 
3.07 
13.13 
3.14 
13.17 
3.23 
13.26 
3.29 
13.29 
3.31 
13.33 
3.35 
13.39 
3.40 
1341 
3.41 
13.43 
3.44 
13.45 
3.46 
13.49 
3.51 
13.52 
3.58 
13.60 
3.62 
13.66 
3.71 
13.77 
3.77 


13.79 


-10.839 
-10.419 
-11.400 
-10.544 
-10.710 
-10.839 
-11.052 
-10.456 
-10.520 
-11.043 
-11.434 
-11.854 
-11.214 
-10.460 
-11.217 
-10.971 
-10.285 
-11.438 
-11.160 
-13.330 
-10.287 
-10.796 
-11.128 
-10.721 
-10.695 
-10.971 
-10.583 
-11.076 
-9.762 
-10.932 
-10.556 
-12.256 
-10.842 
-10.337 
-10.613 
-10.991 
-10.924 
-11.726 
-11.118 
-12.267 
-11.183 
-11.123 
-11.383 
-11.699 


45.895 
47.152 
47.302 
48.150 
47.825 
47.103 
46.779 
50.226 
49.534 
49.257 
46.691 
50.503 
48.980 
47.533 
46.750 
46.645 
47.329 
45.985 
46.964 
44.814 
47.492 
47.329 
47.238 
49.119 
49.949 
49.811 
47.390 
50.365 
47.840 
48.288 
47.167 
44.388 
48.011 
47.014 
47.472 
50.642 
47.714 
45.829 
47.873 
46.602 
48.427 
48.703 
49.396 
47.370 


21.569 
21.254 
22.629 
22.050 
22.214 
22.184 
21.935 
21.636 
21.411 
22.111 
22.734 
22.307 
22.435 
22.030 
22.341 
22.291 
21.771 
22.193 
22.462 
22.339 
21.862 
22.239 
22.107 
22.810 
21.509 
22.131 
22.707 
22.173 
22.218 
21.246 
21.426 
22.138 
21.994 
22.319 
22.323 
23.086 
22.717 
22.652 
21.974 
22.296 
21.819 
21.840 
23.871 
22.375 


46.308 
47.300 
47.853 
47.505 
48.835 
48.071 
47.395 
47.348 
46.744 
46.540 
47.164 
48.101 
47.483 
46.937 
47.376 
47.538 
47.835 
47.106 
47.992 
46.841 
48.343 
48.001 
48.396 
49.195 
47.466 
47.156 
47.160 
48.447 
47.978 
47.097 
48.187 
45.567 
47.156 
47.630 
46.980 
47.650 
48.494 
45.967 
47.478 
47.511 
48.001 
47.743 
47.841 


47.502 


0.34759 


0.05199 
0.13743 
0.18700 
0.03139 
0.02833 
0.13658 
0.04035 
0.29477 
0.03243 
0.05423 
0.10296 
0.09691 
0.04487 
0.04533 
0.11962 
0.09318 
0.01701 
0.17109 
0.04295 
0.13513 
0.01978 
0.02421 
0.00408 
0.06763 
0.04488 
0.33164 
0.12192 
0.05968 
0.04521 
0.19362 
0.07864 
0.05902 
0.05231 
0.15308 
0.00067 
0.00129 
0.00032 
0.00312 
0.00025 
0.00032 
0.00431 
0.00025 
0.00053 


4.478 
4.669 
3.464 
3.948 
4.179 
4.223 
4.274 
4.621 
4.811 
4.069 
3.771 
4.409 
3.819 
4.084 
3.952 
4.207 
4.573 
4.464 
4.222 
4.942 
4.779 
4.319 
4.763 
4.290 
5.116 
4.682 
3.892 
4.674 
4.209 
5.437 
5.503 
4.884 
4.652 
4.433 
4.400 
3.641 
4.195 
3.916 
5.047 
4.813 
5.465 
5.337 
3.114 
4.974 


202102.00005v1 


a 
m 


Xiv 


china 


Chinax ive (ERAT 


4C+55.17 0.896 — 2009/10/31-2009/11/2 13.82 -10.658 47.431 24.289 -9.682 47.897 0.00175 14.94 
GB6 J0654+5042 1.253 2010-1-15 14.01 -10.910 47.030 23221 -10.354 47.586 0.00141 14.77 
PKS 1454-354 1.424 - 14.11 -12.441 45.656 22.441 -11.153 46.945 0.00299 16.00 
PKS 2141+175 0.213 - 15.08 -11.710 44.412 23.918 -11.580 44.542 0.00308 10.00 


2.945 
3.965 
5.290 


4.114 


2 结果 与 分 析 


2.1 同步 与 逆 康 普 峰 的 流量 比 和 多 普 勒 因子 


文 [24] 使 用 均匀 的 单 区 轻 子 模型 进行 分 机 ， 他 们 假设 了 电子 是 以 截断 霸 律 谱 的 形式 注 
入 到 相对 论 性 运动 的 等 离子 体 团 中 ， 注 入 电子 会 随 着 时 间 演 化 。 多 波段 的 辐射 全 部 来 自 这 
个 以 相对 论 性 速度 运动 的 致密 辐射 区 。 他 们 通过 研究 发 现 ， 如果 IC 过 程 是 由 EC 过 程 占 主 
导 ， 则 EC 过 程 贡 献 的 流量 和 同步 峰 流 量 的 关系 近似 如 下 中， 


Up(VFy)Ec 1+а 
—— ző 4 
Vext(VFy) syn ( ) 


HF, Ug = B2/8z 是 人 磁场 能 量 密度 ，Uext 是 外 部 光子 的 能 量 密度 ，6 是 多 普 勒 因子 ，a 是 光 
子 谱 指数 ，B 是 磁场 。 假 设 Uext 是 恒定 的 ""， 常 数 通 常会 离散 数据 ， 并 不 会 影响 相关 性 的 
分 析 。 因 此 ，(2) 式 可 以 预测 IC 峰 流 量 与 同步 峰 流量 的 比值 和 多 普 勒 因子 之 间 的 相关 性 。 
如 果 IC 过 程 是 由 SSC WEES, BBA SSC 过 程 产生 辐射 的 流量 和 同步 辐射 的 流量 之 比 ， 
np DA js 


wF) 2 Es = = 
Tone x 3 (Fr Neo) (M UInXc (ec) X Neo (5) 


其 中 , or 是 汤姆 森 散 射 截 面 ,neo 是 非 热 辐 射电 子 数 密度 的 归 一 化 因子 ,六 是 辐射 区 的 半径 ， 
es 是 同步 峰 分 量 定 义 的 上 限 ，ec 是 高 能 康 普 顿 分 量 定 义 范 围 的 上 限 ， 史 是 转换 后 康 普 顿 - 同 
步 的 对 数 形式 。 此 时 , 我 们 预期 IC 峰 流 量 与 同步 峰 流量 的 比值 和 多 普 勒 因子 之 间 没 有 相关 
性 。 

1 中 展示 了 整个 样本 log(Us (V) gc/ Ue (Vy)syn) 与 log6 的 相关 性 。 斯 皮尔 曼 测 试 显 
示 置 信 水 平 p=0.0003 С Ёр = 0.05) 以 及 相关 系数 R=0.470。 我 们 用 最 小 二 乘法 做 线性 
回归 得 到 的 结果 为 log(UBs(vR,)gc/Uext (VE)syn) = (2.40 + 0.36)1ogó + (—4.76 + 0.43). E 
述 结果 表明 ，FSRQs 中 IC 机 制 主要 由 EC 过 程 主导 。 
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相关 性 分 析 ， 实 线 表 示 最 佳 拟 合 。 


Fig.1 The correlation for the whole sample between log(Us (VE, )c/Vext(vF,) Ба) апа 1028. The solid line is the best linear fitting. 


根据 文 [25] 的 说 法 ， 同 步 峰 峰值 频率 可 能 跟 粒 子 能 被 加 速 到 达 的 最 大 能 量 有 关 。 对 于 


目前 更 为 常见 的 LSP/SP/HSP 分 类 法 E9， 他 们 发 现 统计 上 FSRQs 大 部 分 属于 LSPs， 并 且 


同步 峰 峰 值 频率 集中 在 logvsyn = 13.3Hz， 还 有 逆 康 普 顿 散射 发 生 在 汤姆 进 区 域 (The 


Thomson regime) 。 所 以 在 本 文中 ， 我 们 想 要 确定 是 否 我 们 样本 中 不 同 峰 频 的 FSRQ 也 有 


一 些 物理 上 的 不 同 ， 因 此 我 们 按照 同步 峰 峰 值 频率 集中 在 ~13.3Hz 把 我 们 的 样本 分 成 了 两 
个 子 样本 ， 即 相对 低 同 步 峰 峰 频 种 群 (12.10Hz < logvsyn < 13.05Hz， 包 含 38 个 源 ， 称 低 
峰 频 种 群 ) 与 相对 高 同步 峰 峰 频 种 群 (13.05Hz < logvsyn < 14.11Hz， 包 含 37 个 源 ， 称 高 
峰 频 种 群 ) 。 图 2 中 分 别 展示 了 两 个 种 群 log(Ug F, Јес/ Ue ЈЕ, )ѕуп) 5 logo АНКЕ. 1I 


峰 频 的 种 群 ， 斯 皮尔 曼 测 试 显示 置信 水 平 p=0.013 


以 及 相关 系数 R=0.520。 我 们 用 最 小 二 


we 法 做 线 性 回 归 得 到 的 结 ES A log(Ug(VF,)ec/UVext(VF,) syn) = (2.13 + 0.45)logó t 
(一 4.48 土 0.55)。 高 峰 频 的 种 群 ,斯 皮尔 曼 测 试 显示 置信 水 平 p=0.008 以 及 相关 系数 R=0.433。 


我 们 用 最 小 二 乘法 做 线性 回归 得 到 的 结果 为 


log(Ug(vF,)gc/Uexi (УЕ, ) syn) = (2.72 + 


0.58 1028 + (—5.09 土 0.69)。 上 述 分 析 结 果 表 明 ， 在 分 成 两 个 子 样本 进行 研究 得 到 的 结 


与 整个 样本 的 结果 没有 明显 差异 ， 表 明 在 FSRQs B 


ЕЖ, IC 机 制 主 要 由 EC 过 程 主 导 。 
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图 2 两 个 种 群 og(Us(vE,)gc/Uext(VR,)syn) 与 log6 的 相关 性 分 析 。 左 图 为 低 峰 频 种 群 (12.10 На Лови у, <13.05 Hz) 包含 38 + 
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源 ， 右 图 为 高 峰 频 种 群 (12.10 Hz<logveyn<13.05 Hz) 包含 37 个 源 。 实 线 表 示 最 佳 拟 合 。 
Fig.2 The correlation for two populations between log(Us VE,)Ec/ Ue (Rss) and logó. The left panel is the low frequency peak 
populations (12.10 Hz<logvsyn «13.05 Hz) containing 38 sources, and the right panel is the high frequency peak populations (13.05 


Hz<logvsyn «14.11 Hz) containing of 37 sources. The solid line is the best linear fitting. 


22 种 子 因子 


在 上 一 节 的 研究 中 , 我 们 发 现 FSRQs 的 IC 辐射 都 是 由 EC 过 程 主导 的 , 在 EC 过 程 中 
存在 多 种 软 光 子 起源 , 例如 吸 积 盘 (accretion disk, AD) , 宽 线 区 Cbroad-line region, BLR) 
和 人 尘埃 环 (dust torus, DT) 等 。 为 了 进一步 调查 其 软 光子 的 来 源 。 我 们 使 用 文 [17]Harvey 
提出 的 自 举 法 “bootstrapping) 去 计算 低 、 高 峰 频 种 群 的 观测 的 种 子 因子 分 布 seed factor 
distribution) (BULA 1) 。 计 算 种 子 因 子 的 表达 式 为 "1: 


1/2 


SF = 10810( SE log10 (s. 22 x 104 —— Gauss, (6) 


VIC,22 


其 中 ，Uext 是 喷 流 外 部 光子 能 量 密度 (在 CGS 制 下 ) ; eo 是 外 部 光子 的 特征 能 量 〈 以 静止 
怨 子 的 能 量 为 单位 ) ; CD 是 康 普 顿 主导 (CD = (vLv)icpeak/ (VLV)synpeak， 即 IC 峰 的 峰 
值 光度 与 同步 峰 的 峰值 光度 的 比值 ,以 10 为 单位 );vsyni3 是 同步 峰 的 峰值 频率 (以 1013 为 单 
BED ; viczz 是 IC 峰 的 峰值 频率 〈 以 1022 为 单位 ) ВІВ 和 DT 理论 上 期 望 的 种 子 因子 分 
1 ， 是 通过 (6) 式 的 外 部 光子 场 及 外 部 光子 特征 能 量 进行 估算 的 。 我 们 采用 Harvey 给 出 
的 范围 ， 即 : SFgig = 3.29 + 0.11; ЗЕрт = 3.92 + 0.11. 

我 们 使 用 双 侧 柯 尔 莫 蕊 洛 夫 检验 (two-sided Kolmogorov-Smirnov test) 对 各 个 样本 的 
种 子 因子 分 布 进行 了 正 态 性 检验 。 我 们 发 现在 低 、 高 峰 频 种 群 中 均 不 能 拒绝 正 态 性 假设 。 
在 低 峰 频 种 群 中 ， 拒 绝 正 态 性 假设 的 显著 性 为 0.376 (р=0.71) ; 而 在 高 峰 频 种 群 ， 拒 绝 正 
态 性 假设 的 显著 性 为 0.016 (р=0.99) 均 小 于 26 的 置信 水 平 。 因 此 ， 这 些 样本 的 种 子 因 子 
分 布 服从 正 态 分 布 。 与 Harvey 的 研究 一 致 。 

基于 自 举 法 ， 我 们 从 低 、 高 峰 频 种 群 中 观测 的 种 子 因子 分 布 分 别 构 建 了 观测 种 子 因子 
的 大 样本 〔 低 峰 频 种 群 38X10 个 ， 高 峰 频 种 群 37X10 个 ) ,分 别 从 大 样本 中 构建 10* + 
包含 相同 源 数目 的 奉 代 样 本 。 使 用 替代 样本 ， 分 别 估算 了 不 同 种 群 中 观测 的 种 子 因子 分 布 
的 中 位 数 CObsqeaian? 2. 26 置信 区 间 ， 即 : 

低 峰 频 种 群 ，Obsmedian = 4.20+0:25 (7) 
高 峰 频 种 群 ，Obsmeqian = 4.43*021 (8) 

图 3 中 展示 了 两 个 种 群 的 种 子 因子 分 布 直方 图 ， 并 且 标注 了 理论 上 期 望 的 宽 线 区 种 子 

因子 和 人 尘埃 环 种 子 因 子 分 布 。 从 图 3 上 ， 可 以 发 现 直方 图 峰值 处 靠近 尘埃 环 的 位 置 。 
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—— Dust toru SF 
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Histogram Frequency 


Seed factor 


频 种 群 (13.05 Hz<logy, 


图 3 两 个 种 群 的 种 子 因子 分 布 直方 图 。 左 图 是 低 峰 频 种 群 (12.10 Hz<logvsyn<13.05 


Seed factor 


Hz) ， 包 含 38 个 种 子 因子 ， 右 图 是 高 峰 


——- Broad-line region SF 
Broad-line region SF range 
一 一 Dust toru SF 
~ Dust toru SF range 
WE Histogram Frequency 


yn<14.11 Hz) ， 包 含 37 个 种 子 因子 。 直 方 图 用 黑色 区 域 展示 。 蓝 色 实 线 和 虚线 分 别 表示 理论 上 的 宽 线 
区 种 子 因子 及 范围 。 黄 色 实 线 和 虚线 分 别 表 示 理 论 上 人 尘埃 环 的 种 子 因子 及 范围 。 直 方 图 的 区 间 是 由 
自动 选择 的 。 从 两 幅 直 方 图 上 ， 可 以 发 现在 直方 图 的 峰值 处 靠近 尘埃 环 的 位 置 。 


matplotlib. pyplot.hist 


Fig.3 Histogram of two populations observed seed factor. The left panel is the low frequency peak populations (12.10 


Hz<logvsyn «13.05 Hz) containing the seed factor of 38 sources, and the right panel is the high frequency peak populations (13.05 


Hz«logv,,, «14.11 Hz) containing the seed factor of 37 sources. The histogram is plotted in black. The expected broad-line region seed 


factor and range are shown in blue, with solid and dashed lines, respectively. The expected dust torus seed factor and range are shown 


in yellow, with solid and dashed lines, respectively. The histogram was binned using the auto option for the matplotlib. pyplot. hist. 


As can be seen from this histogram, the seed factor distribution peaks close to the seed factor distribution of dust torus. 


种 子 因子 分 布 的 直方 图 ， 没 有 包含 估计 的 不 确定 性 。 为 了 包含 估计 的 不 确定 性 ， 基 于 
自 举 法 实现 了 低 、 高 峰 频 种 群 的 核 密度 估计 Ckernel density estimation, KDE) 。 运 用 


Silverman’s Rule?” КЮЕ 中 的 带宽 自动 进行 了 平滑 处 理 中 
度 估 计 。 在 图 4 中 ， 


表明 软 光 子 很 可 能 
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图 4 两 个 种 群 的 核 密度 估计 。 左 图 是 低 峰 频 种 群 (12.10 Hz<logvsyn<13.05 Hz) , 
(13.05 Hz<logvsyn<14.11 Hz) ， 包 含 37 RFI. 5 


的 种 子 因子 分 布 不 相交 ， 表 明 软 光子 可 能 
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seef factor 
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SAR EP BER MER, 


seef factor 
含 38 个 源 的 种 子 因 子 ， 右 图 是 高 峰 频 种 群 
灰色 区 域 是 为 了 视觉 辅助 。 蓝 色 实 线 和 虚 


线 分 别 表示 期 望 上 种 子 因子 及 范围 ; 黄色 实 线 和 虚线 分 别 表示 期 望 上 的 尘埃 环 种 子 因子 及 范围 ; 绿色 实 线 和 虚线 分 别 表 示 观 测 


上 的 种 子 因 子 中 位 数 及 26 置信 区 间 。 外 观 上 , 核 密度 估计 近似 一 个 正 态 分 布 。 在 20 的 置信 水 平 下 ， 
示 均 不 能 拒绝 正 态 性 假设 ， 从 左 到 右 拒绝 正 态 性 假设 的 显著 性 分 别 为 0.376 (р=0.71) 和 


К 


双 侧 柯 尔 莫 戈 治 夫 检 验 显 


0.016 (р=0.99) 。 


Fig.4 Kernel density estimate of two populations. The left panel is the low frequency peak populations (12.10 Hz<logvsyn <13.05 Hz) 


containing the seed factor of 38 sources, and the right panel is the high frequency peak populations (13.05 Hz<logvsyn «14.11 Hz) 


containing the seed factor of 37 sources. The black solid line is the kernel density estimation, and the gray area is the visual aid. The 


expected broad-line region seed factor and range are shown in blue, with solid and dashed lines, respectively. The expected dust torus 


seed factor and range are shown in yellow, with solid and dashed lines, respectively. The median of observed seed factor and 26 
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confidence interval are shown in green, with solid and dashed lines, respectively. The KDE has the general appearance of a normal 
distribution. A two-sided Kolmogorov-Smirnov test indicates that normality cannot be rejected significantly (using a rejection 


significance threshold of 26), with rejection significance of 0.376 (p = 0.71) and 0.016 (p = 0.99) from left to right, respectively. 


为 了 进一步 检验 观测 的 种 子 因子 分 布 是 与 宽 线 区 还 是 尘埃 环 位 置 对 应 ， 我 们 分 别 计算 
了 低 、 高 峰 频 种 群 中 观测 种 子 因 子 分 布 与 理论 上 期 望 的 种 子 因子 分 布 之 间 的 拒绝 中 位 数 相 
等 的 显著 性 。 在 标准 差 方面 ， 对 于 低 峰 频 种 群 ， 拒 绝 的 显著 性 为 : 


CCObsmedian — BLRmedian) = 8.280 (9) 

o(Obsmedian 一 DT median) = 2.550 (10) 
而 在 高 峰 频 种 群 ， 拒 绝 的 显著 性 为 ; 

O(ObSmedian ^ BLRmedian) = 10.390 (11) 

o(Obsmedian 一 DT median) = 4.660 (12) 


因此 ， 低 峰 频 中 观测 的 种 子 因子 分 布 与 理论 上 期 望 的 DT 种 子 因子 分 布 的 拒绝 显著 性 
最 小 ， 表 明 低 峰 频 种 群 中 软 光子 的 起 源 应 该 很 可 能 是 由 尘埃 环 主导 的 。 而 高 峰 频 种 群 中 虽 
然 观测 种 子 因子 分 布 和 DT 种子 因子 分 布 的 显著 性 比 与 BLR 种 子 因子 分 布 的 拒绝 显著 性 小 ， 
但 是 KDE 显示 出 明显 差异 ， 表 明 软 光子 的 起 源 对 于 人 尘埃 环 或 许 不 是 单一 的 位 置 ， 软 光子 
可 能 有 BLR 和 DT 的 贡献 。 


3 讨论 与 总 结 


我 们 研究 发 现 无 论 是 整个 样本 还 是 按 数 量 均 分 的 两 个 子 样 本 ， 在 
log(UBp(vR,)Ec/UBLR(VE,)syn) 和 log6 之 间 的 相关 性 方面 ， 图 1 与 图 2 没有 表现 出 明显 的 差 
异 。 这 就 表明 ， 轻 子 模型 的 框架 下 ，FSRQs 中 的 IC 机 制 很 有 可 能 由 EC 过 程 主导 。 图 1 
中 包含 了 75 个 FSRQs， 而 图 2 左 和 右 分 别 包 含 了 38 和 37 个 FSRQs， 图 上 的 数据 分 布 两 
者 均 存 在 一 定 的 弥散 ,这 从 某 种 程度 上 会 削弱 我 们 的 结论 。 随 着 越 来 越 多 观测 设备 的 加 入 ， 
=== 未 来 可 以 收集 到 更 多 的 样本 进一步 的 验证 轻 子 模型 框架 下 ， 高 能 辐射 的 主导 过 程 。 

种 子 因子 的 计算 基于 多 波段 SEDs 的 拟 合 ， 要 求 宽带 辐射 没有 明显 的 光 变 现象 "”。 我 
们 使 用 Harvey 提出 运用 自 举 法 估算 观测 种 子 因子 分 布 的 方法 去 计算 我 们 样本 的 种 子 因 子 。 
图 3 的 直方 图 展示 了 观测 的 种 子 因子 分 布 ， 该 直方 图 没有 包含 估计 的 不 确定 性 ， 图 4 的 核 
密度 估计 包含 了 估计 的 不 确定 性 .从 图 3 可 以 发 现 ,观测 的 种 子 因子 分 布 集中 在 峰 频 位 置 ， 
而 图 4 的 KDE 曲线 同样 展示 了 相同 的 结果 ， 并 且 展 示 了 由 我 们 的 样本 计算 得 到 观测 的 种 
子 因子 分 布 与 Harvey 理论 上 期 望 的 种 子 因 子 分 布 之 间 的 差异 。 从 图 4 (AE) 可 以 发 现 ， 观 
测 的 种 子 因子 分 布 在 统计 上 与 DT 种 子 因 子 分 布 的 范围 相交 ， 表 明 在 低 峰 频 的 种 群 中 软 光 
子 很 可 能 由 尘埃 环 主导 。 从 图 4〈 右 ) 可 以 发 现 ， 观 测 的 种 子 因 子 分 布 超过 了 理论 上 期 望 
的 宽 线 区 和 尘埃 环 种 子 因子 分 布 的 范围 ， 表 明 在 高 峰 频 种 群 中 DT 很 可 能 不 是 主导 软 光子 
的 单一 位 置 ， 可 能 有 BER 和 DT 共同 的 贡献 。Harvey 的 研究 发 现 FSRQs 的 辐射 区 与 DT 
的 位 置 在 统计 上 一 致 ， 表 明 软 光子 主要 是 DT 的 贡献 。 与 Harvey 的 研究 结果 相 比 ， 我 们 的 


样本 在 不 同 种 群 中 显示 出 了 不 同 的 结果 。 图 3 的 直方 图 分 别 包含 38 个 和 37 个 种 子 因子 ， 


因此 统计 上 计算 的 种 子 因 子 仿 少 ， 会 削弱 我 们 的 结论 ， 和 希望 未 来 进一步 增加 同时 性 样本 ， 


并 且 使 用 不 同方 法 获得 的 观测 量 来 作 进一步 分 析 验 证 
本 文 从 统计 上 研究 了 75 个 平 谱 射 电 类 星体 中 软 光 子 的 起 源 。 通 过 相关 性 分 析 , 结果 表 


о 


明 平 谱 射 电 类 星体 的 IC 机 制 很 可 能 由 EC 过 程 主 导 。 运 用 种 子 因子 分 布 进一步 检验 了 EC 
过 程 中 软 光 子 的 来 源 ， 结 果 表 明 不 同 的 种 群 在 统计 上 有 不 同 的 结果 。 
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Abstract: To verify dominating processes in the inverse Compton mechanism and the 
origin of soft photons in the dominating processes. We collected multi-wavelength 
(quasi-)simultaneity data of 75 flat-spectrum radio quasars from the literature. We fitted the 
spectral energy distribution using logarithmic parabolic functions, then analyzed three 
correlations for the whole sample, the subsamples of two bins, respectively. The research results 
are summarized as follows. (1) We find significant correlations of 
log(Ug(VF, ec/Uext(VF,)syn) and logó for the whole sample and subsamples of two bins, 
respectively. Our results show that the inverse Compton mechanism of the flat-spectrum radio 
quasar is preferentially dominant by the external Compton process under the leptonic model. (2) 
Based on the Bootstrap to calculate observed seed factor distribution, we apply this method to 
constrain the origin of soft photons of the dominant process of inverse Compton mechanism. The 
Observed seed factor distribution suggests that soft photons are likely to be dominated by dust 
torus for populations with relatively low synchrotron peak frequencies. But, we cannot determine 
the origin of soft photons of the dominant process of inverse Compton mechanism for 
populations with relatively high synchrotron peak frequencies. 

Key words: Jet-radiation mechanisms; Active-galaxies; Flat-spectrum radio quasar; Inverse 


Compton scattering 


